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• Relativistic corrections communicate the binding energy of a bound state to its kineti c mass. This mec hanism 

is reviewed and used to explain anomalous results of Collins, Edwards, Heller, and Sloan ( |hcp-lat/9512026| ), which 
compared rest and kinetic masses of heavy-light mesons and quarkonia. 



1. INTRODUCTION 

Last year Collins, Edwards, Heller, and 
Sloan studied heavy Wilson quarks || 
with the improved action of Shcikholcslami and 
Wohlert 0. In lattice units the heavy-quark 
mass rriQa typically exceeded unity, a regime in 
which the numerical results require a nonrelativis- 
tic interpretation [|||] just as in NRQCD {]|J. 
Ref. presents a test of the nonrelativistic in- 
terpretation that removes kinematic effects and 
focuses on a dynamical effect — the binding en- 
ergy. The results of this test were unexpectedly 
anomalous. The aim of this paper is to explain 
why and to offer a remedy. 

Let us begin with some basics to define nota- 
tion. As a function of momentum p the energy of 
a state X can be written 



E x (p) = M 
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where the rest mass is defined by M\ 
the kinetic mass is defined by 



M 2 



d 2 E 
d PiJ P =o 



(1) 

E(0), and 



(2) 



Below the states can be quarks, Q and q, and 
mesons QQ, Qq, and qq. Usually Q is assumed 
heavier than q. 

On a relativistic mass shell 
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(3) 



In this paper the lower-case mx denotes the ex- 
act, physical mass, whereas upper-case Mix de- 
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note the result of a (possibly approximate) cal- 
culation. In the mass-dependent renormalization 
of ref. H it is possible to adjust the action's cou- 
plings to recover eq. (||), to a specified accuracy. 
With the Wilson or Sheikholeslami- Wohlert ac- 
tions, however, perturbation theory shows that 
Miq ^ M 2 q (except as rriQa —►()). In nonrel- 
ativistic systems this is acceptable, provided one 
adjusts the bare mass until the kinetic mass takes 
its physical value ^|||, just as in ref. |||7]]. 

The quark state makes sense at most in pertur- 
bation theory. In a nonperturbative world, one 
would like to carry out the tuning with a bound 
state, e.g. a meson X = Qq whose masses M^^, 
M 2 Q q , etc, can be computed with Monte Carlo 
methods. Let us define the binding energy B by 



M lQq = M 1 Q + M lq + B 1 , 
M 2Qq = M 2 Q + M 2g + B 2 . 



(4) 



To make a precise definition of the binding ener- 
gies, one requires a precise definition of the quark 
masses in eq. (|^). Here it is enough to take the 
rest and kinetic masses of the free theory (<?q = 0) 
with the same bare quark masses. Computing 
bound-state splittings through rest masses makes 
sense only if B\ is the experimental binding en- 
ergy, to sufficient accuracy. Similarly, tuning the 
meson mass to the kinetic mass makes sense only 
if also B 2 is the experimental binding energy, to 
sufficient accuracy. 

As shown below, the anomalous result of ref. ||] 
is sensitive to B 2 — B\ . There are two keys to un- 
derstanding it. First, one must be careful about 
the qualifying phrase "to sufficient accuracy" in 
the preceding paragraph. Second, one must to 
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know how field theories communicate the bind- 
ing energy to a bound state's kinetic mass. 

Sect. H recalls the diagnostic test of ref. (!]]. 
Sect. | assesses the cutoff effects of the binding 
energy in light-light and heavy-light mesons, and 
in quarkonium. The mechanism for generating 
the "kinetic binding energy" is reviewed for a 
relativistic (continuum) gauge theory in sect. 
and generalized to lattice gauge theory in sect. pT 
Sect. | draws a few conclusions. 

2. THE TEST 

Let us abbreviate SM := M 2 - Mi and SB := 
B 2 - B x . Ref. § introduces 



2SM Qq - (SM Q 



QQ 



SM, 



(5) 



Comparison with eqs. (Q) shows that the quark 
masses drop out, leaving 



I = 



25B r 



(SB, 



QQ 



SB, 



(G) 



If the lattice action(s) of the quarks were suffi- 
ciently accurate, all SBs, and hence /, would van- 
ish. (/ vanishes trivially when tuq = m q , even if 

The numerical results of ref. Q are shown in 
fig. [I]. The "inconsistency" / is negative, and |/| 
tends to increase with increasing tuq . To explain 
both the sign and the magnitude, below I shall 
derive an expression for SB. 

3. CUTOFF EFFECTS ON SB 

Before presenting the analytical result for SB, 
it is useful to anticipate the order of magnitude 
of SB in each meson — light-light, heavy-light, and 
quarkonium. On this basis it turns out that the 
quarkonium SBqq dominates the numerator in 
eq. (1). 

3.1. Light-light SB qq 

The binding energy is O(Aqcd)- With the 
Sheikholeslami-Wohlert action, B\ and B 2 both 
suffer from lattice artifacts of 0(a™aAQCD)- 
With the tree-level improvement of used by 
ref. pj, n = 1. (With the Wilson action n = 0.) 




Figure 1. Plot of the binding-energy "inconsis- 
tency" / vs. the heavy-light meson's kinetic mass 
aM 2 Q q , for vtlq 3> m q . Adapted from ref. Q. 

There is no reason for the artifacts to be identi- 
cal, so SBq q is 0(a™aAQ CD ). This is numerically 
small, so SBqq can be neglected below. 

3.2. Heavy-light SB Qq 

The binding energy is again O(Aqcd)- The 
light quark suffers lattice artifacts as above, but, 
when rriQ ^ 1 — as in fig. |l|, the heavy quark also 
suffers from (smaller) effects of 0(Aq CD o/mQ). 
Again, even though there is no reason for artifacts 
in Bi and B 2 to cancel, one sees that SBq q is 
numerically negligible. 

3.3. Quarkonium SBqq 

The binding energy is now 0(mQV 2 ), where v 
denotes the relative Q-Q velocity. In this nonrela- 
tivistic system, the velocity is a pertinent estima- 
tor of cutoff effects @,|[. The rest mass is 0(t>°), 
so the action would need absolute accuracy of 
0(v 2 ) to obtain relative accuracy of 0(v 2 ) in B\. 
Both the Wilson and Sheikholeslami-Wohlert ac- 
tions achieve this. On the other hand, the kinetic 
mass multiplies an 0(v 2 ) effect, so the action 
would now need an absolute accuracy of 0(v ) to 
obtain relative accuracy of 0(v 2 ) in B 2 . Neither 
the Wilson nor the Sheikholeslami-Wohlert action 
achieves this H ; with either of them, one can only 
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hope for 0(v ) relative accuracy in B 2 . The er- 
ror in B 2 , and hence in SBqq ^ 0, is 0(mQV 2 ), 
which is significantly larger than the previous two 
estimates. 

4. BREIT EQUATION 

For nonrelativistic systems the binding-energy 
discrepancy can be worked out quantitatively, fol- 
lowing a textbook nonrelativistic expansion of 
QED |1. This section verifies in a relativisti- 
cally invariant theory that B 2 = B\ = B. The 
next section then turns to the lattice theories, 
which break relativistic invariance. For conve- 
nience, these two section assume that even the 
"light" quark q is nonrelativistic. 

At leading order the quark-anti-quark inter- 
action arises from one-gluon exchange diagrams. 
Evaluating these diagrams and developing the 
nonrelativistic expansion, one obtains a Hamilto- 
nian H = ttlq + m q + H 2 + H& for the quark- 
anti-quark system. The leading nonrelativistic 
dynamics are given by 



Ho = - 



P 



Q 



2rrin 
™Qq 



2m q 
2(i 



-V(r) 
V(r), 



(7) 



where V(r) = —CFCt/r; r — xq — x q , P and p 
are center-of mass coordinates and momentum; 
fx = (rriQ 1 + m^ 1 )^ 1 is the reduced mass; and 



M Q q = m Q- 
are 



The first relativistic corrections 



(P|) 2 
8m% 



(P 2 q ? 
8m q 



V 2 (r;p Q ,p q ), 



(8) 



where the non-local potential V 2 is given by 
Breit's equation M. It takes the form 



V 2 (r;pQ,p q ) = const x <5^(r) 



V(r) 



1 - 



PQ-P q + r nrjPQiPqj 



(9) 



spin-dependent . 



The spin-dependent terms and the terms propor- 
tional to 5^ 3 \r) are not important here. Full 



details are given in §§ 83-84 of ref. [|S|. To- 
gether with the (p 2 ) 2 terms in H4, the exhib- 
ited part of V 2 is responsible for modifying the 
bound-state kinetic mass from Mq q — ttiq + m q 
to Wq 9 = mQ + m q + B (as required by Lorentz 
invariance). 

To proceed one must re-write H4 in center-of- 
mass momenta and collect terms quadratic in the 
total bound-state momentum P. In the bound 
state, combinations of the internal momentum p 
and relative coordinate r can be replaced by ex- 
pectation values. Collecting all terms, the bound- 
state kinetic energy becomes 



2m 



Qq 



1 - 



{T + V) 



2M Qq V M Qq 
P P 



(10) 



where T — p 2 j2\i is the internal kinetic energy. 
By the virial theorem the second line vanishes. 
Thus, to consistent order in p/M the leading rela- 
tivistic corrections H4 generates the right binding 
energy B 2 = (T + V) =: B for the bound-state 
kinetic mass. 

More generally, higher-order relativistic effects 
trickle down to bound-state properties as follows: 
the correction of 0(v e ) provides the 0(v e ~ k ) con- 
tribution to bound-state properties of 0(v k ). 

5. LATTICE GENERALIZATION 

On a hypercubic lattice there can be two cor- 
rections to the kinetic energy 

(P 2 ) 2 



E{p) 



8M| 6 



1 W40 3 p\ 



for each of p = pq , p q . Here 

( v 1/3 . 



dp 2 dp 2 



p=0 



and 



W4 := 



1 d 4 E 



4 dpi 



p=0 



4Mf 



(11) 



(12) 



(13) 



Unless the action has been improved further than 
the Shcikholcslami-Wohlert action, M4 ^ M 2 and 
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w 4 ^ 0, cf. Appendix A of ref. ||. These lattice 
artifacts filter through to B 2 — just as above — 
through the terms proportional to PiPj{piPj). 

The spatial gluon generates the spin- 
independent contribution proportional to V{r) in 
eq. ([|); on the lattice the nuance is that the ki- 
netic mass appears in the bracket. On the other 
hand, the temporal gluon generates more compli- 
cated terms, but they either depend on spin or 
are proportional to 5^ (r). So to work out an ex- 
pression for i? 2 , it is enough to maintain eq. (|^), 
but with masses M 2 Q q and /i 2 built from M 2 q 
and M 2q . 

The calculation of the binding energy differ- 
ence 5B follows the steps leading to eq. One 
finds an expression that is too cumbersome to 
present here. In an S wave, however, it can be 
simplified because (piPj) — \5ij(p 2 ). Then 




/'2 



M h 



Mi, 



(14) 



+ Aa 2 f i 2 {M 2 Q w 4Q + M 2 q w iq ) 



This is the main new result of this paper. Note 
that, as one would have anticipated, the expres- 
sion vanishes when w^x = and M^x — M 2 x- 

With an estimate of (T) from potential mod- 
els H and the lattice masses of the right-most 
point in fig. [l], I find 



SB 



QQ 



2M, 



-0.5. 



(15) 



The agreement with the Monte Carlo results is 
surprisingly good. 

6. CONCLUSIONS 

The origin of the anomaly observed in ref. Q| 
is the usage of an action accurate only to 0(v 2 ). 
Thus the relative error in the binding energy B 2 
of the bound-state kinetic mass is of order 
ttiqv 2 /mQV 2 = 1. Meanwhile, the usual binding 
energy B\ is indeed valid to leading order in v 2 . 
The test quantity I cleverly isolates B 2 —B\ 1 and 
thus exposes an inconsistency of 0(1). 



By examining how (approximately) relativistic 
field theories generate B 2 , this paper explains the 
results found last year |l|. Moreover, the analy- 
sis makes the remedy plain: the anomaly is not 
expected to appear if quarkonium properties are 
computed with an action improved through 0(u 4 ) 
(or higher). In particular, one requires Mi = M 2 
and u>4 = 0. 

Most published applications of ref. use a suf- 
ficiently accurate action 0. Ref. (To) even re- 
marks that 0(i; 4 ) accuracy is essential for a con- 
sistent determination of the 6-quark mass from 
the T spectrum. For four-component fermions 
the details required for 0(v 4 ) accuracy in quarko- 
nium have appeared more recently J5[. 
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